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ABSTRACT 

The ASCA and BeppoSAX spectra of the Circinus Galaxy and NGC 1068 are anal- 
ysed and compared with photoionization models based on cloudy. In the case of 
Circinus, a single, mildly ionized reflector can account for the line spectrum, while 
in NGC 1068 at least three different reflectors (with different ionization states) are 
needed. We suggest that the reflector in Circinus, and the low ionized one in NGC 1068 
are the inner and visible part of the material responsible for the X-ray absorption. 
With this assumption, we estimate for the inner radius of the absorber a value of 0.2 
pc for Circinus and of a few parsecs for NGC 1068. 

Key words: Line: formation - galaxies: individual: Circinus Galaxy, NGC 1068 - 
galaxies: active - X-rays: galaxies 



1 INTRODUCTION 

The heavy obscuration of the nuclear X-ray emission in 
Seyfert 2 galaxies allows to observe and study components 
that in type 1 sources are diluted into invisibility by the nu- 
clear radiation. Among these components, particularly inter- 
esting are those due to reflection of the nuclear X-rays from 
circumnuclear matter, as they permit in principle to map 
the reflecting regions and probe their physical and chemical 
properties through the study of the continuum and espe- 
cially of the emission line spectra (e.g. Netzer 1996; Matt, 
Brandt & Fabian 1996; Krolik 1999 and references therein). 

If the absorbing matter is Compton-thick (i.e. Nh > 
a' 1 = 1.5 x 10 24 cm" 2 , see Matt et al. 2000 for the discus- 
sion of the general properties of Compton-thick sources), the 
nuclear radiation is obscured up to at least ~10 keV, and all 
astrophysically important X-ray lines (including iron lines) 
emitted by the reflecting matter can be studied. 

In this paper we discuss the ASCA and BeppoSAX X- 
ray spectra of the two brightest and best studied Compton- 
thick sources, i.e. the Circinus Galaxy and NGC 1068. In 
the case of NGC 1068, we present for the first time results 
from the A04 100 ks ASCA observation (results published 
so far by Ueno et al. 1994, Iwasawa, Fabian & Matt 1997, 
and Netzer & Turner 1997, were obtained with the PV 40 
ks observation). Our aim is to derive the properties of the 
reflectors in these sources. To achieve this goal, we have 



extensively used the photoionization code cloudy* (version 
90.05; Ferland 1997), to construct theoretical models to be 
compared with the observations. 

The paper is organized as follows: Sec. 2 summarizes 
previous X-ray observations of the two sources; Sec. 3 de- 
scribes the data reduction, while Sec. 4 reports the spectral 
analysis and Sec. 5 compare the results with the theoreti- 
cal models. In each section we will discuss first the Circinus 
Galaxy and then NGC 1068. Finally, the results will be sum- 
marized in Sec. 6. 



2 PREVIOUS X-RAY OBSERVATIONS 
2.1 The Circinus Galaxy 

The Circinus Galaxy, one of the closest AGN and X-ray 
brightest Seyfert 2 galaxies, was observed for the first time 
in X-rays during the ROSAT All Sky Survey (Brinkmann, 
Siebert & Boiler 1994). Later on, it was observed by ASCA; 
the spectrum is dominated by a pure Compton reflection 
component, with a very prominent iron Ka line and sev- 
eral other lines from lighter elements (Matt et al. 1996). 
The lines from the lighter elements have energies indicating 
significant ionization of the emitting material (Matt et al. 
1996; Netzer et al. 1998; Sako et al. 2000). The BeppoSAX 
observation, besides confirming the ASCA findings below 

* http://www.pa.uky.edu/" gary/cloudy/ 
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10 keV, detected the nuclear emission, seen through very 
thick (4xl0 24 cm" 2 ) absorbing matter (Matt et al. 1999; 
Guainazzi et al. 1999). 

2.2 NGC 1068 

The history of the X-ray observations of NGC 1068 is longer 
and richer than that of the Circinus Galaxy. The source was 
observed by Einstein (Monier & Halpern 1987) and EX- 
OSAT (Elvis & Lawrence 1988); the lack of any measured 
intrinsic absorption was interpreted by those authors as evi- 
dence for reflected emission of an obscured nucleus. This sce- 
nario was confirmed by the Ginga detection of a strong iron 
line (Koyama et al. 1989). At softer energies, the ROSAT 
HRI discovery of a very extended (~13 kpc) emission indi- 
cated that a large part of the radiation in this band is related 
to the host galaxy rather than the nucleus. 

ASCA (Ueno et al. 1994; Iwasawa, Fabian & Matt 
1997; Netzer & Turner 1997) confirmed both the thermal- 
like emission in soft X-rays and the reflection component 
at higher energies. The very rich line spectrum and the 
complex iron line (including a neutral or low ionized line 
as well as He- and H-like lines) indicate that more than 
one reflector is present. Smith, Done & Pounds (1993) pro- 
posed to identify the cold reflector with the inner surface 
of the putative molecular torus, an idea further pursued by 
Matt, Brandt & Fabian (1996) and Iwasawa, Fabian & Matt 
(1997). BeppoSAX confirmed the complexity of the reflec- 
tion in NGC 1068 by separating at least two reflecting con- 
tinua, one cold and one ionized (Matt et al. 1997; Guainazzi 
et al. 1999). The lack of any transmitted component implies 
that the absorbing matter should have a column density ex- 
ceeding 10 25 cm" 2 . 

Finallly, evidence for flux and spectral variability be- 
tween the two BeppoSAX observations is discussed by 
Guainazzi et al. (2000). 



3 DATA REDUCTION 
3.1 The Circinus Galaxy 

In the spectral analysis discussed below we will mainly use 
ASCA data (observation performed on 1995, August 14), 
but we will also take advantage of the BeppoSAX results 
(observation performed on 1998, March 24). 

Details on the BeppoSAX observation and related data 
reduction can be found in Matt et al. (1999) and Guainazzi 
et al. (1999). Regarding ASCA, we have verified the previous 
results (Matt et al. 1996) using the Revision2 data process- 
ing and the latest calibration. It should be noted that the 
SIS0 spectrum appears to be affected by a gain problem: 
e.g., the iron K line at 6.4 keV is shifted to a higher en- 
ergy (see below for details) and a similar energy shift is seen 
for the instrumental Ni Ka line at 7.5 keV. We have checked 
the instrumental Ni Ka line at 7.5 keV for the two SIS data 
taken from a source-free region on the detector field of view 
and found that the SIS0 data show an energy shift similar 
to that seen in the source spectrum while the SIS1 data 
show the line at the right energy. We therefore believe that 
the line energies derived from the SIS1 are more reliable. 
Also, a recent Chandra HETG observation of the Circinus 



galaxy found a narrow line at 6.4 keV (R. Sambruna, private 
communication) . 



3.2 NGC 1068 

In the spectral analysis discussed below we have used both 
ASCA (observation performed on 1996, August 15) and Bep- 
poSAX (observations performed on 1996, December 30 and 
1998, January 11) data. 

Details on the BeppoSAX observations and related data 
reduction can be found in Matt et al. (1997) and Guainazzi 
et al. (1999, 2000). 

The ASCA data presented here are from a long ob- 
servation carried out in the AO-4. The data reduction was 
performed in the standard procedure using FTOOLS ver- 
sion 4.2 and the latest calibration available. The exposure 
time for both SIS detectors is 118 ks, respectively. 



4 SPECTRAL ANALYSIS 
4.1 The Circinus Galaxy 

We have re-analysed the ASCA SIS1 spectra of the Circinus 
Galaxy to search for emission lines, taking advantage of the 
results on the continuum obtained with BeppoSAX (Matt 
et al. 1999). We have fitted the spectrum with the following 
model: 

F — [Rc + PLs + emissionlines] exp{— Up^A^sr} (1) 

where Rc is the cold reflection component (model pexrav; 
the power law photon index of the nuclear illuminating con- 
tinuum has been fixed to the value of 1.56 obtained with 
BeppoSAX, Matt et al. 1999), PLs is a power law compo- 
nent describing the soft X-ray emission, and Nh is the cold 
absorption along the line of sight, which includes absorp- 
tion from our own Galaxy plus possible contributions from 
the Circinus Galaxy itself. The primary, heavily absorbed 
continuum discovered by BeppoSAX has not been included 
in the model because it does not affect the ASCA energy 
band. We have added as many emission lines as required 
by the data, inspecting by eye the residuals after any sin- 
gle addition. All lines have been assumed to be narrow for 
simplicity, even if, as discussed in Sec. 5.1, some of them are 
very likely a blend of different lines with comparable fluxes. 
With all lines included, the reduced \ is 0.8 for 141 d.o.f.. 
The best fit parameters for the continuum are listed in Ta- 
ble 1, while Fig. 1 and 2 show the best fit models for ASCA 
and BeppoSAX. 

We will focus our interest on the line spectrum in Sec. 
5.1: the best fit line parameters are listed in Table 3, along 
with the most likely identifications. SIS1 data are systemat- 
ically different from those of SIS0: the line energies (even if 
always at least marginally consistent with the other instru- 
ment) are all shifted to higher values, while the equivalent 
widths (except for the first one) are all much lower than 
those of SIS0. The reason for this discrepancy is not clear. 
The EWs measured by SIS0 may be more reliable, since 
this instrument is considered to be best calibrated at lower 
energies, while the correct line energies are probably those 
measured by SIS1 (see discussion in Sec. 3.1). The value of 
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Figure 1. ASCA SIS1 data and best model for Circinus. 
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Figure 2. BeppoSAX LECS+MECS data and best fit model for 
Circinus. 



the reduced x 2 for the fit of the continuum is similar for the 
two instruments. 

Upper limits on the equivalent widths of the ionized 
iron lines are 30 eV for the 6.7 keV Fexxv line, and 340 eV 
for the 6.96 Fexxvi line. 

The photon index of the soft power law is not very well 
constrained, i.e. T s = 2.661q4 3 , while the column density of 
the absorbing matter is 6.6(±2.1) x 10 21 cm~ 2 . (The fit to the 
BeppoSAX data, after fixing the same value for the column 
density, gives a flatter power law with T s = l.l^oll)- A 
thermal plasma model for the soft excess does not provide 
a better fit to the data. The temperature is 1.84to. 2 ? keV, 
while the abundance is constrained to be less than 0.7 times 
the solar value. As the thermal model cannot account for the 
observed lines, we do not discuss it any further. The 0.5-2 
keV luminosity of the soft component is around 7xl0 39 erg 
cnf 2 s _1 . The IR/soft-X-ray ratio is similar to that of the 
sources in the sample of powerful infrared galaxies studied 
by Iwasawa (1999). 



NGC1068 (ASCA) 




Energy (keV) 

Figure 3. ASCA SIS0+SIS1 data and best fit model for NGC 
1068. 



4.2 NGC 1068 

As already noted by Iwasawa, Fabian & Matt (1997) and 
Guainazzi et al. (1999), the X-ray spectrum of NGC 1068 
is much more complex than that of the Circinus Galaxy. 
The continuum is composed by at least three components: 
a thermal-like one, probably associated with the extended 
starburst observed by ROSAT (Wilson et al. 1992) and dom- 
inating the soft X-ray band; and at least two reflectors, one 
cold and one ionized (Iwasawa, Fabian & Matt 1997, Matt 
et al. 1997). Superimposed on these components are many 
emission lines, most of them associated with the ionized re- 
flector^), even if the neutral Ka line is clearly visible. 

We adopted for the continuum the same model of 
Guainazzi et al. (1999): a single temperature thermal plasma 
model (mekal in Xspec), a cold reflection component 
(pexrav) and an ionized reflector component, assumed to 
be a simple power law with the same spectral index of the 
power law illuminating the cold matter, which in turn has 
been fixed to 2.13, as measured by BeppoSAX (Guainazzi et 
al. 1999). The column density of the absorbing matter has 
been fixed to the Galactic value of 3.5 xlO 20 cm -2 (Dickey 
& Lockman 1990). Then, as in Circinus, we added emission 
lines until required by the data. The best fit parameters for 
the continuum are given in Table 2, giving a reduced % 2 °f 
~ 1.1 for 582 d.o.f.. They are consistent with those obtained 
by Guainazzi et al. (1999) from the BeppoSAX data. 

The free parameters in the fit, apart from the energies 
and the normalizations of the lines, are the temperature and 
the abundance of the mekal model: T = 0.7 ± 0.1 keV and 
A z = (l.3±g;|) x 10" 1 . 

A fit to the BeppoSAX data (Fig. 4) with the same 
model gives a similar value for the reduced x 2 , but the 
mekal parameters are not consistent with the ASCA ones: 
T = 0.40±o.oi keV and A z = (l-6±g;|) x 10~ 2 . Two more 
lines are needed for the BeppoSAX data with respect to 
those found with ASCA, at energies of 0.55 and 8.2 keV. 
This is the result of the larger band coverage of the instru- 
ments aboard BeppoSAX. The lines at 0.55 e 8.2 keV have 
been both observed and discussed by Guainazzi et al. (1999). 
The identification of the second one with the K/3 transition 
of Fexxvi at 8.2 keV (Guainazzi et al. 1999) allows to fix 
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Table 1. Best fit parameters for the continuum of Circinus. The values with the symbol * were kept fixed during the fit. 





ASCA 


BeppoSAX 




powerlaw 


mekal 


powerlaw 


mekal 


reduced x 2 


0.8 (141 d.o.f.) 


0.8 (140 d.o.f.) 


1 (92 d.o.f.) 


1 (91 d.o.f.) 


r 


1.56* 


1.56* 


1.56* 


1.56* 


cut-off (keV) 


56* 


56* 


56* 


56* 


N H (cm' 2 ) 


(6.6+21) x 1Q 21 


(5.7tH) x 10 21 


6.6 x 10 21 * 


6.6 x 10 21 * 


T s 


2 6fi+ - 62 
Z - DO -0.43 




i 7^+0.28 
1 -'°-0.92 




kT (keV) 




1 84+ 1 ' 2 




349+11 


A z 




(i-s±5:i) x 10-1 




o q + G.52 



NGC1068 (BeppoSAX) 




Table 2. Best fit parameters for the continuum of NGC 1068. 
The values with the symbol * were kept fixed during the fit. 



Energy (keV) 



ASCA BeppoSAX 


reduced x 


~1.1 (582 d.o.f.) 


~1.3 (111 d.o.f.) 


r 


2.13* 


2.13* 


T s 


2.13* 


2.13* 


N H (cm' 2 ) 


3.5 x 10 20 * 


3.5 x 10 20 * 


kT (keV) 


0.7 ±0.1 


0.40+;;- 


A Z 


(l.3t° ;i) x lO-i 


(l.6+° ; 3 2 ) x IO" 2 


redshift 


3.793 x 10~ 3 * 


3.793 x IO" 3 * 



Figure 4. BeppoSAX LECS+MECS+PDS data and best fit 
model for NGC 1068. 



its energy and calculate the ratio against the Ka of Fe xxvi, 
which is 0.17^+°. 
The line at 0.55+ 

CK02 keV confirms one of the most in- 
teresting results of Guainazzi et al. (1999), i.e. the detection 
for the first time of the O VII line. It is worth remarking that 
this line is in excess of what expected from the thermal-like 
component, and therefore it is likely to be produced in one 
of the reflectors. This line was not observed in previous ob- 
servations because of the low energy resolution (ROSAT) or 
the narrow band coverage (ASCA). The upper limit set by 
the BBXRT observation suggested that the abundance of 
the oxygen could be less than the cosmic one (Marshall et 
al. 1993; Netzer & Turner 1997). Instead, the clear presence 
of this line in the BeppoSAX data overcomes this difficulty 
and, as we will see in Sec. 5.2, puts strong constrains to 
the number and the characteristics of the reflectors in NGC 
1068. 

All emission lines are listed in Tables 6 and 7: we will 
discuss them in detail in Sec. 5.2. 



5 COMPARISON WITH MODELS 
5.1 The Circinus Galaxy 

5.1.1 The model 

The unambiguous presence of a Compton reflection compo- 
nent in the ASCA (Matt et al. 1996) and BeppoSAX (Matt 
et al. 1999) spectra indicate the presence of thick circum- 
nuclear matter. From the ratio between the reflected and 
direct components, it is likely that the covering factor of 
this matter is fairly large (Matt et al. 1999, Guainazzi et 
al. 1999, Matt et al. 2000). It is then natural to assume 
that the reflecting and absorbing media are one and the 
same. We have therefore calculated with cloudy the re- 
flected spectrum from a X-ray illuminated material with a 
column density of 4.3xl0 24 cm -2 (Matt et al. 1999), adopt- 
ing the dimensionless ionization parameter U, defined as: 



U = 



Q 



4:nr 2 cn 



(2) 



where Q is the total number of ionizing photons per unit 
time reaching the surface of a cloud at a distance r from the 
source, while n is the hydrogen number density. With the 
nuclear X-ray spectrum and luminosity of Circinus (Matt 
et al. 1999), we obtain (cgs units): 



U ■ 



1.04 x 10 40 



(3) 



© 0000 RAS, MNRAS 000, 000-000 



The circumnuclear X-ray reflectors in NGC 1068 and the Circinus Galaxy 5 



Prom the measurements of the dynamical mass, Risaliti et 
al. (1999) have set an upper limit of ~20 pc on the external 
radius of the cold matter, 7?, in agreement with recent esti- 
mates from IR observations, which give a value of R ~ 16 pc 
(Ruiz et al. 2000). Thus, assuming r <C R, the density is 
~ 8.7 x 10 4 cm' 3 , and the inner surface of the cloud is 
r ~0.11C/~2 p C . As described below, we have found a solu- 
tion with C/~0.25, which implies a value for r of ~0.2 pc, 
roughly consistent with the water maser measurements of 
Greenhill et al. (1997) and with the value of 0.5 pc esti- 
mated from the IR observations whose R we adopted in our 
calculations. Then the assumption r <C R is, a posteriori, 
justified. (It must be noted that the above determination of 
r is based on the simplyfing assumption of constant density 
along the torus. A negative density gradient would, however, 
be probably more realistic. If this is the case, this value of r 
must be regarded as an upper limit.) 

We searched for the simplest possible solution, i.e. a 
single reflector. The iron line undoubtedly comes from low 
ionized atoms, while lighter elements are fairly ionized. We 
therefore searched for the highest value of U consistent with 
iron less ionized than xvn, to ensure that the iron line is 
still at ~6.4 keV, and that the resonant trapping (e.g. Matt, 
Fabian & Ross 1993) is not effective. We found such a so- 
lution with U =0.25. The ionization structures for iron, sul- 
phur, silicon and magnesium are shown in Fig. 5. 

As a consistency check, we refitted the data with an 
ionized reflector (i.e. pexriv in Xspec) instead of the cold 
one. The ionization parameter used by pexriv has a differ- 
ent definition: 

f = — erg cms' 1 (4) 
nr 1 

where L is the 5 eV — 20keV luminosity, assuming the 
continuum as a simple power law (Magdziarz & Zdziarski, 
1995). The conversion factor between the two definitions is 
^ — 249 erg cm s _1 . Thus, the ionization parameter used in 
our model corresponds to £ = 249 • 0.25 ~ 62 erg cms' 1 . 

We first checked if reflection from ionized material could 
explain the soft excess, and then we fitted ASCA and Bep- 
poSAX data without the additional soft power law. In both 
cases, the fit is unacceptable, giving a reduced y 2 greater 
than 2. We then fitted the data including the power law 
component, and fixing £ = 62. The ionization parameter has 
been fixed because it was not well constrained: the value ob- 
tained with cloudy is more reliable, as pexriv calculates 
the ionization equilibria in a simplified way (Done et al., 
1992). In these conditions, the ASCA fit has a reduced y 2 
almost identical to the one obtained with pexrav. The soft 
power law photon index is V s = 2.18 + . ) 'gg, while the ab- 
sorbing column density N H = (3.6±i'|) x 10 21 cm~ 2 . The 
flux in the power law is now a factor 2 lower than for the 
cold reflector. The results for BeppoSAX data are roughly 
the same, but again with a photon index (T s = 0.7lio'_o?) 
which is not consistent with the ASCA one. 

Finally, we have estimated the solid angle of the 
reflecting matter which is observable by comparing the 
observed F a = 3.1(±0.3) x 10~ 4 ph cm' 2 s' 1 ~ 3.2 x 
\Q' 12 erg cm' 2 s' 1 ) , and calculated iron line fluxes. We find 
a value of 0.2x47r. 



Table 3. X-ray emission lines in Circinus from ASCA data (SIS1 
and SISO). The EWs are calculated against the reflection con- 
tinuum (pexriv). All line profiles have been assumed to be 6- 
functions. 



SIS1 SISO 



E(keV) 


EW(eV) 


E(keV) 


EW(eV) 


Id. 


n oo+0-05 
at "-0.05 


Q1n +910 


n Q9+0-04 
U - y -0.04 


nnnn + 2500 
2500 -1500 


Ne ix 


i o 4 +0.05 


645+ 180 


-, , s +0.04 
1,oo -0.03 


™ + 310 
oou -310 


Mgxi 


1 - oo -0.02 


910-22 


1 - 8 '-0.02 


' Ju -190 


Si xm 


2 16+ - 03 


580±|?8 


2 40+ 06 
z -^ u -0.04 


•jnc+130 


Si xm 


9 -9+0.11 
z - ,z -0.32 


ioq+130 


d - UU -0.30 


g +80 
uu -60 


S xv/xvi 


o qc+0.07 
J - uo -0.07 


170+ 110 

±lu -110 


o 99 + 0.08 
° -0.07 


1 "50+ 90 
iou _100 


Ar xvn 


6 41+ 002 


2130l 2 ^ 


6 49+0-02 
°- 4la -o.oi 


1D(U_ 160 


Fe Kol 


7 00+ - 04 

' - uu -0.04 


640 + 170 


6 98 +0 ' 1() 


9 60 +l00 
ZDU -100 


Fe K/3 



5.1.2 The line spectrum 

The spectrum is rich in emission lines, and it is likely that in 
most cases the observed lines are actually a blend of two or 
more lines with comparable intensities. In Tab. 4 we report 
the observed equivalent width of the lines (measured with 
respect to the reflected continuum only) compared with the 
total EW of the most important lines likely contributing to 
any single blend, as calculated by cloudy. The agreement 
between observed and calculated one is good for the lines 
with energies above 1.85 keV provided that the abundance 
is significantly higher than the solar one, i.e. Az ~ 3 for the 
SISO data (A z ~ 4 for the less reliable SIS1 data) (However, 
it should be noted that Oliva et al. (1999) found, instead, 
evidence for lower than solar metal abundances). 

This value for Az is still consistent with the Fe Ka line, 
which does not depend linearly with the abundance, but 
only logarithmically (Matt, Brandt & Fabian 1996), and it 
is strongly dependent on the other elements' abundance. For 
Az 2± 3, the EW of the iron line increases by only a factor of 
~ 1.4, in agreement with the observations. Since there are 
not detailed calculations of this effect for the other lines, we 
assume, like cloudy apparently does, that their dependence 
on the abundance is almost linear; in any case, the derived 
Az should be considered as a lower limit. It should be also 
noted that the model predict two other weak line systems, 
centered at at 2.09 and 2.56 keV, which are not observed, 
not surprisingly given the very low equivalent widths. 

In the table we did not include the line at ~ 7keV, 
which could be identified either as the K/3 from neutral iron 
(7.06 keV), or the line from Fexxvi (6.97 keV). Matt et al. 
(1996) proposed a blend of these two lines, because they 
observed a very large EW. This kind of solution would be 
unacceptable for our model, because it would imply the pres- 
ence of a highly ionized reflector for the Fe xxvi line to be 
produced. If we consider the SIS1 data, we face the same 
problem: the EW of the line, if it has to be identified with the 
Kf3 only, should be ~ | the Ka, while the observed one is 
much more intense, 0.3±0.1. If, instead, we analyse the data 
from SISO, the ratio Ka/K/3 becomes consistent with the 
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Figure 5. Ionization structures in Circinus for Fe, S, Si and Mg. 



expectations, being 0.15 ±0.06. On the other hand, cloudy 
does not seem to include the Fe Kf3 line in its calculations, 
so we could not compare it with the data. 

The comparison of the lines observed in the ASCA spec- 
trum at 0.82 and 1.34 fceV with the cloudy model is made 
difficult by the large number of lines present, in the calcula- 
tions, at those energies. These lines are mainly K lines from 
Ne and Mg and L lines from Fe, but contributions from Ca, 
Mn and other elements are also present. When folded with 
the ASCA and BeppoSAX response matrices, these lines 
are blended to a pseudo-continuum, which contributes for 
about 10% of the observed soft excess. A detailed compar- 
ison between models and data for these lines must await 



high resolution spectra, as those provided by the gratings 
onboard Chandra and XMM-Newton. 

Given the very simple assumptions, the agreement be- 
tween observations and model is remarkable. At the resolu- 
tion of a CCD, therefore, a single, mildly ionized material 
seems able to explain the Circinus X-ray spectrum. 

Finally, as a consistency check, we have also calculated 
the optical emission lines from such a reflector, and com- 
pared with those observed by Oliva, Marconi & Moorwood 
(1999) in the central, unresolved region of the source. As 
there may be other regions contributing to the optical spec- 
trum besides the X-ray reflector, we are contented with not 
overproducing the observed lines. This goal is achieved as- 
suming a dust-to-gas ratio typical of the interstellar matter, 
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Table 4. X-ray emission lines produced in our model for Circinus. 
The energy in the first column is the mean value (weighted on the 
flux) of all the lines that contribute to that blend. The last column 
shows the ratio between the observed EW (SISO data) and that 
calculated by cloudy. All the lines are K transitions. The Ka 
lines which dominate each blend are in boldface. 



£(keV) 


EW(eV) 


Id. 


E(keV) 


EW a 
EW C 


1.85 


245 


Si xiii 

Mgxn 
Mgxn 


1.85 

1.87 
1.91 


3.2 ± 0.8 


o no 


1 n 
1U 


Si xiv 
Si xm 


2.17 




2.37 


52 


Si xm 
Si xm 

Sxv 


2.29 
2.38 
2.45 


5.9 ±2.5 


2.56 


5 


Si xiv 
Si xiv 
S xvi 


2.50 
2.60 
2.63 




2.99 


15 


Sxv 
Sxv 
Ar xvii 


2.88 
3.02 
3.08 


40 +5.5 


6.45 


1300 


FeKa 


6.45 


1.3 ±0.1 



and provided that the temperature of the reflecting mat- 
ter is larger than ~10 J K, otherwise the [O IIIJA5007 line 
is overproduced (an agreement with this line is obtained 
for T ~2.5xl0 5 K). The calculated flux of the Hj3 is in- 
stead a few times lower than observed, as for other intense 
lines (Tab. 5). It is worth noting that dust is expected to 
be present, since the reflecting surface of the torus is well 
beyond the evaporating radius, which is (Barvainis 1987): 

r ev ~L% 2 pc~ 0.07 pc (5) 

if we assume Ltot — Lib, — 5 x 10 43 erg s' 1 (Siebenmorgen 
et al. 1997). It must be recalled that dust may survive in an 
ionized medium provided that the electron temperature is 
lower than 10 6 K, see Draine & Salpeter (1979). The pres- 
ence of dust, however, is in disagreement with the very low 
Av/Nh ratio for this source. A possibility is that the size 
of the grains is significantly larger than in the interstellar 
medium (Salvati & Maiolino 2000; Maiolino et al. 2000); an- 
other possibility is that the optical/IR and X-ray absorbers 
are physically distinct, which is well possible as the optical 
and X-ray emission regions are certainly different. 

5.2 NGC 1068 

5.2.1 The model 

Differently from the Circinus Galaxy, a single reflector can- 
not explain the X-ray spectrum of NGC 1068. In fact, it 
is impossible to have from a single region lines from neu- 
tral and He-like/H-like iron atoms. Indeed, Guainazzi et al. 
(1999) showed that at least three reflectors are needed. We 
used cloudy to derive the physical properties of these re- 
flectors, which in the following will be called the cold, warm 
and hot reflectors. 



Table 5. Optical emission lines in Circinus (Hf3 = 100). 



A (A) 


ESO NTT 


CLOUDY 


Id. 


6563 


298 


70 


Ha 


5007 


1025 


260 


Oin 


4959 


320 


90 


Om 


4861 


100 


26 


Hf) 


4686 


41 


50 


He ii 


3869 


136 


10 


Ne in 



The need for a third reflector arises from the presence 
of the O VII fluorescent line. It can be produced in the cold, 
thick reflector (which, as in Circinus, can be naturally iden- 
tified as the torus itself), but constraining the ionization pa- 
rameter to be too low to allow the production of all the other 
lines, except for the 'neutral' Fe Ka. On the other hand, in 
the hot reflector, responsible for the He- and H-like iron 
lines, all the lighter elements are fully stripped. The Mg, Si 
and S lines must therefore be produced in a third reflector, 
the warm one. 

Let us investigate the characteristics of the three re- 
flectors in more detail. The illuminating radiation has been 
modeled as a power law with photon index 2.13 (Matt et 
al. 1997). The cold reflector is made of low ionized material 
and can be identified with the inner walls of the torus. We 
thus used a column density of 10 26 cm -2 (Matt et al. 1997) 
and an electron density of n e — 10 6 cm' 3 , corresponding 
to an outer radius of 30 pc. The best ionization parameter 
is U — 0.1, with a temperature of 10 K . The ionization 
structures for some elements are shown in Fig. 6. 

The warm reflector produces most of the remaining ob- 
served lines, except for the Fexxv and Fexxvi lines. The 
best ionization parameter is loglJ = 0.7, and the tem- 
perature is 10 K. The ionization structures for some el- 
ements are shown in Fig. 7: the oxygen is fully stripped, 
while Fexxv and Fexxvi have low abundances. We used 
a column density of 10 21 cm -2 and an electron density of 
n e — 10 6 cm -3 (see below). 

The lines from the ionized iron come from the hot re- 
flector, with a column density similar to that of the warm 
one but an ionization parameter of logU = 2.1. The ion- 
ization structure for some elements are shown in Fig. 8: all 
elements but iron are fully stripped. 

The values adopted for the column densities of the warm 
and the hot reflectors have two justifications. The first con- 
strain, which concerns the hot reflector only, is from the EWs 
of the ionized iron lines which, as we will see in the next sec- 
tion, are in good agreement with the results of Matt, Brandt 
& Fabian (1996) with N H ~ 10 21 cm' 2 . The other constrain, 
which applies to the sum of the two reflectors, comes from 
the amount of the ionized reflection continuum, considering 
its dependence from the nuclear luminosity of the source: 

L r ~ L n ■ tc ■ f < L n ■ tc (6) 

where / = ^ < 1 is the covering factor and tc = 
1 ^" n is the Compton optical depth. Assuming the ob- 
served ASCA value for the warm reflection component, 
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Figure 6. Ionization structures in the cold reflector of NGC 1068 for Fe, O, Si and Mg. 



L r — 0.6 ■ 10 41 erg s 1 , we can have an estimate of the nu- 
clear X-ray luminosity: 

r ^ 9 ' 1043 -1 

L n > 21 erg s (7) 

Since the X-ray luminosity of NGC 1068 is estimated to be 
around 10 44 erg s' 1 (Iwasawa, Fabian & Matt 1997) (corre- 
sponding to a bolometric luminosity of a few 10 45 erg s^ 1 
adopting the SED of Elvis et al. (1994), roughly in agree- 
ment with the estimate of Pier et al. 1994), we see from (7) 
that we should use Nh > 10 21 cm~ 2 . This is also the value 
adopted by (Guainazzi et al. 1999). 

Finally, from the analysis of the EW of the lines it can 
be shown (see next section) that the continua produced in 



the warm and in the hot reflectors should have column den- 
sities of the same order of magnitude, which justifies the 
above assumptions. 

We can estimate the inner radius of the torus in NGC 
1068 similarly to what have been done for Circinus. The 
difficulty here is that the nuclear luminosity and the density 
of the matter are not known. With the best guess values for 
these two quantities we have: 

r ~ 3.9 LH 2 44 n- 1/2 pc (8) 

For a discussion of the probable X-ray luminosity see Iwa- 
sawa et al. (1997). For the density, we can only derive a 
lower limit. In fact: 
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Figure 7. Ionization structures in the warm reflector of NGC 1068 for Fe, S, Si and Mg. 



n e = — j^- > —j— ~ 3 x 10 4 • (-Rio) 1 cm 3 (9) 
R R 

with the outer radius of the torus R in units of 10 pc. The 
(9) is a lower limit both because the column density could 
be grater than 10 25 cm -2 , and because an upper limit on 
R of ~ 100 pc (Risaliti et al. 1999) can be derived from 
the dynamical mass (Thatte et al. 1997) and the mass of 
the gas itself measured from CO observations (Heifer 1997). 
Our estimate for the inner radius is in agreement with the 
water maser measurements of Gallimore et al. (1996), which 
have shown the presence of a keplerian disk with an inner 
radius r in > 1.3pc and outer radius r out ~ 2.5pc. 

Following the same considerations made for Circinus, 
we reanalysed ASCA and BeppoSAX data using the pexriv 



model instead of pexrav. The fits are of the same sta- 
tistical quality as the previous ones, with no significant 
changes in the values of the other parameters. The ioniza- 
tion parameter which minimizes the \ 2 m BeppoSAX data 
is £ = 28ljg erg cm s _1 , marginally consistent with the one 
calculated by cloudy (£ ~ 5 erg cm s~ : the conversion fac- 
tor is different from the one we used in Circinus because the 
power law slope is different). 

Finally, we can estimate the visible fraction of the re- 
flecting surface in NGC 1068, comparing the fluxes of the 
Fe Ka in our model with the observed one. We obtain: 

0~O.l (10) 

using r ~ 3.9 pc. 
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Figure 8. Ionization structures in the hot reflector of NGC 1068 for Fe and S. 



Table 6. Fe K complex of NGC 1068 from ASCA data. The values 
with the symbol * were kept fixed during the fit. The EWs are 
calculated against different continua, as explained in Sec. 5.2.2. 



E (keV) 


EW (eV) 


Ident. 


6.40* 


900+ 130 
auu -l00 


Fe Ka 'neutral' 


6.70* 


43501™° 


Fe xxv 


6.96* 


Q40+940 

y4U _82o 


Fe xxvi 


7.06* 


100* 


Fe Kj3 'neutral' 



5.2.2 The line spectrum 

Let us first discuss the iron line complex. In Tab. 6 the Fe 
K lines from ASCA data are shown. All line energies were 
kept fixed in the fit, as well as the Fe Kj5 flux (g of Ka). 
The EWs are calculated against the reflection continuum of 
the material that produces the line: Fe Ka and K/3 against 
the cold reflection component (pexra(i)v model in Xspec, 
cold reflector in cloudy), while Fexxv and Fexxvi against 
the ionized component (power law in Xspec, hot reflector 
in cloudy). 

We compared the observational results for the iron lines 
with the work of Matt, Brandt & Fabian (1996) (from now 
on MBF96), because it is specifically devoted to this prob- 
lem, while we continue using cloudy to calculate the rela- 
tive abundances of the iron ions. 

As discussed by MBF96, the Fexxv line is emitted both 
by resonant scattering and recombination; in the moderately 
thick case (for the resonant line) of Nh = 10 21 cm -2 , its EW 
can be expressed as: 

EWfcxxv ^ (2100 • fxxv -A z ) + (1000- fxxv ■ A z ) eV(ll) 
(see Fig. 5 of MBF96) where fxxv = 0.57 is the ion abun- 



dance calculated by cloudy. We must remember that this 
equivalent width has to be multiplied by a factor X < 1 (see 
previous section) to be compared with the observed one: 



A z -X 



4350 



2.5 



(12) 



(3100 ■ 0.57) 

The EW of the Fe xxvi line can be expressed in a similar 



way: 

EWfeXXVI 



(U00-fxxvi-A z )+(^0-fxxvi-A z ) eV(13) 



Using fxxvi = 0.22, as calculated by cloudy, and the 
estimate of Az ■ X (eq.12), we have: 



EW F exxvi = 990 eV 



(14) 



in good agreement with the observations. The value in (12) 
is, therefore, consistent with both the ionized iron lines pro- 
duced in the hot reflector. 

In Tab. 7 all the other lines observed by ASCA are 
shown. Almost all these lines are produced in the warm re- 
flector: therefore, in the table are also reported the EWs 
multiplied by a factor 1 — X = 0.25 which is our estimate 
for the contribution of the warm reflector alone (see below). 
From (12), this implies Az — 3.3. This is the value used for 
the metal abundances in our model: in this case, being the 
material optically thin, the dependence of the EWs from Az 
is linear. Like in Circinus, emission lines have been blended 
when necessary: the results are summarized in Tab. 8. 

Our choice for X is justified by the comparison of Tab. 
8 with Tab. 7. The choice of X = 0.75 and A z = 3.3 is 
the one which gives a good agreeemnt between model and 
observations. With these parameters, the comparison is sat- 
isfying apart from a deficit in the Si lines. Regarding the 
lines at energies lower than 1.85 keV, the model spectrum is 
too rich to allow for a useful comparison with the data, as 
already noted for Circinus. 

The O vil line is produced, as explained in the previous 
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Table 7. X-ray emission lines in NGC 1068 from ASC'A data 
(SIS1 and SISO). The values obtained by dividing for by the 
factor 1 — X = 0.25 (see text) are also shown. 



iii (Kev ) 


EW (eV) 


Ident. 




observed 


1 - X = 0.25 




1 07+ ' 2 


190+H 




Fe l 


i-20±g;l 


910+210 




FeL 


i.35t° .; 


290+^0 


11601-° 


Mg xi 


l.8llS;« 


330+« 


1320l 1 1 6 2°o 


Si xm 


1 96+ ' 3 
i.yo_ 2 






Si xiv 


2.38i;^ 


290 -50 


1160+ 240 

llou -200 


S xv - Si xiv 


2.57+° j 




740+ 320 


S xvi 


s.io±8:g 


1501™ 


ouu_ 240 


Arxvn 



Table 8. X-ray emission lines produced in our model for NGC 
1068. The energy in the first column is the mean value (weighted 
over the flux) of the lines that contribute to that blend. The lines 
with the symbol * are contributes from the hot reflector. All the 
lines are K transitions. The Ka lines which dominate each blend 
are in boldface. 



"E(keV) 


EW (eV) 


Id. 


E (keV) 






Si xm 


1.85 


1.85 


1070 


Mgxn 
Mgxn 


1.87 
1.91 


2.01 


640 + 960* 


Si xiv 


2.01 


2.29 


860 


Si xm 


2.29 






S xv 


2.45 


2.50 


1210 


Si xiv 
Si xiv 


2.50 
2.55 


2.60* 


1000 


Si xiv 


2.60 


3.11 


570 


Ar xvii 
Sxv 


3.08 

3.14 



section, by the cold reflector: its observed EW, calculated 
against the pexriv continuum is 4000l 2 ggg eV. This value is 
in good agreement with our model, which predicts 4000 eV, 
if we assume, as it seems reasonable, Az = 3.3 also for the 
cold reflector. Incidentally, it should be noted again that 
the assumption of linear dependence of EWs from Az in 
optically thick materials is not necessarily a good one. 

The Fe Ka at 6.4 keV is also produced in the cold reflec- 
tor. The equivalent width measured by ASCA (900+. eV) 
is rather low when compared with the results of MBF96, 
taking also into accont the metal overabundance discussed 
above, which implies a further increase of the line flux by 
a factor of ~ 1.4. This correction contributes to make the 
EW of the 'neutral' Fe Ka not consistent with ASCA data: 
using a mean inclination angle, our model predicts a value 
of ~ 1600 ■ 1.4 ~ 2200 eV, a factor two larger than the ob- 
served one. On the other hand, BeppoSAX data are in agree- 



ment with ASCA, giving an equivalent width just slightly 
higher (~ 1100 eV). This problem has no obvious solutions 
at present. 



6 SUMMARY 

The results of our analysis may be summarized as follows: 

For the Circinus Galaxy, a single, mildly ionized and 
optically thick reflector can explain the line spectrum and 
most of the continuum. Only the soft excess remains un- 
explained, and it may be related to the starburst region. 
Recently, Sako et al. (2000) have modeled the line spectrum 
of Circinus with the sum of a cold and thick plus an ionized 
and thin reflectors. At the spectral resolution of ASCA, both 
models are tenable. We only note that our model is simpler, 
making use of only one, instead of two, reflectors. High reso- 
lution observations with the gratings onboard Chandra and 
'KMM-Newton will certainly clarify this issue. 

In NGC 1068, instead, at least three reflectors are 
needed, one low, one mildly and one highly ionized. (A 
thermal-like component, probably associated with the ex- 
tended emission discovered by ROSAT/HRI, see Wilson et 
al. 1992, dominates in the soft X-ray band). In compari- 
son with the Netzer & Turner (1997) model, our scenario 
for the line spectrum is necessarily more complex, because 
the detection of the O vn line by BeppoSAX requires one 
more reflector. On the other hand, contrary to Netzer & 
Turner(1997) and following MBF96 and Iwasawa, Fabian & 
Matt (1997), we identify the low ionized reflector with the 
very thick absorber (Matt et al. 1997), so reducing by one 
the number of regions needed to explain the overall spec- 
trum. 

We identify the reflector in Circinus, and the low ion- 
ized one in NGC 1068, with the inner and visible part of 
the material responsible for the X-ray absorption. With 
this assumption, we estimate for the inner radius of the ab- 
sorber a value of 0.2 pc for Circinus and of a few parsecs for 
NGC 1068. 
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